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Abstract 

We have considered the BCS-BEC crossover at T = in the frame- 
work of a lattice fermion model (the extended Hubbard model). We have 
analyzed the case of on-site hard core repulsion {U — oo) and intersite 
attraction {W < 0) in a 3D system with extended s-wave pairing sym- 
metry. The mechanisms of density-driven and attraction-strength-driven 
crossover have been analyzed. We have obtained the ground state phase 
diagram. The analyzed quantities are: chemical potential, order parame- 
ter, coherence length and fraction of condensed particles. 

1 Introduction 

High- Tc superconductivity has been widely discussed since its discovery by Bed- 
norz and Miiller in 1986 [1]. A microscopic mechanism which could explain 
the characteristic properties of High- Tc superconductors has not been found 
yet. However, there has been broad agreement that the physics of the BCS- 
BEC crossover is crucial for High- Tc superconductivity [2]. The effects of the 
crossover are clearly seen in the behaviour of the energy gap and the chemical 
potential. The so-called pseudogap has been observed in the excitation spectra 
of many unconventional superconductors [3] . 

In 2004 Jin's group experimentally observed a condensation of pairs of 
fermionic atoms in the region of the BCS-BEC crossover [4]. A trapped gas 
of fermionic atoms of ''"K was cooled to suitably low temperatures of the order 
of 5 • 10^^ K and thereafter the interatomic interaction strength was controlled 
via the Feshbach resonance. 

In this paper we analyze the physics of the BCS-BEC crossover in High- Tc 
superconductors by means of a lattice fermion (the extended Hubbard) model 
for the (3D) simple cubic lattice system with hard core on-site repulsion and 
intersite attraction. This paper consists of three parts. The first part gives a 
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discussion of the extended Hubbard model in the Hartree-Fock approximation. 
The second part presents the properties of selected quantities and the numerical 
results. In sec. Ill we summarize the discussion. 



2 Model 

The model Hamiltonian is the extended Hubbard model with on-site (U) and 
intersite (W) interactions: [5], [6], [7]: 

where rijo- — cl^Cic, tij - hopping integral. Having transformed Hamiltonian (1) 
to the reciprocal lattice and using the Hartree-Fock approximation, the following 

system of equations for the superconducting order parameter (A/j = Ao+7/cA^) 
in the case of the s* pairing symmetry is obtained; 



1 + UMT) UMT) \fAo 



0, (2) 



where: ^i(T) = j^^g F,{T), MT) = ^ l,F,iT), cj,^{T) ^ j^Y., llF,{T), 
Fq = {2Eq)~^tanh{f3Eq/2), ~ 2{coskxax + coskyUy + cosfczCz), - lattice 
constant in the j-th direction (we set = 1 in further considera- 

tions), z - coordination number. The equation for the chemical potential takes 
the form: 

n-l = -^Y.e-kFk{n (3) 
fe 

where ek = Ck — fi, n - electron concentration. Solving the above equations 
simultaneously, one can find the dependence of the energy gap and the chemical 
potential on the electron concentration. Taking into account only the hopping 
between the nearest neighbours, the electron dispersion is: 

€k = —2t{coskx + cosky + cosfcz). (4) 

The evolution from the weak coupling limit (BCS) to the limit of the strong 
coupling (BEC) takes place when we decrease the electron concentration or 
increase the interaction. According to the Leggett criterion [8], the BCS-BEC 
crossover takes place when the modified chemical potential: 

A = /u - n{U/2 + Wo) (5) 

drops below the lower band edge. 

When [/ — > 00, the system of equations for the order parameter takes the 
form (provided that Aq 7^ 0, A^ 7^ 0): 

$i(T)Ao + $2(T)A^ = 0, (6) 
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Figure 1 : The ground state phase diagrams {U = oo) for the 3D case (left) and 
the 2D case (right). 



_I5^$2(T)Ao+fl-^$^(T)) A^ = 0. (7) 



z \ z 



Because of the fact that we consider the case of U/t = oo, the on-site pairing 
is forbidden. One can show it by calculating the anomalous corelation function 
[9]. 



3 Results 

We have performed an analysis of the BCS-BEC crossover for the model with 

hard core repulsion. The analyzed quantities arc: chemical potential, order 
parameter, coherence length and fraction of condensed particles. 

Fig. 1 shows the ground state (T = 0) phase diagram for the 3D case (left) 
and the 2D case (right; this case has been analyzed in detail in Rcf [9]). Both 
in the 2D and the 3D case there exists a critical interaction value for which 
a bound state forms in the empty lattice (when n ^ 0). For the 2D case 
\Wc\/t = 2. At the same time, it is the minimal interaction value for which a 
transition to the normal state takes place. The very existence of the transition 
to the normal state at T = at suitably low values of the intersite interaction is 
interesting - it indicates that quantum fluctuations arc very strong. For the 3D 
case the critical interaction value for which a bound state forms in the n — * 
limit is |Wc|/f » 5.874 but the critical interaction value for which a transition 
to the normal state takes place is |W^c|/i = 2. Therefore, as oposcd to the 2D 
case, the minimal interaction value for the emergence of superconductivity is 
different than the critical interaction value for the emergence of a bound state 
in the empty lattice. In both cases the range of occurrence of the local pairs 
phase widens when we decrease the electron concentration or we increase the 
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Figure 2: (left) Order parameter at the Fermi level vs. electron concentration; 
the arrows show the BCS-BEC crossover points, (right) The dependence of 
the modified chemical potential on the electron concentration and the intersite 
interaction |W| (inset); the horizontal dashed line denotes the band bottom. 



interaction. Because of the fact that we consider the case oiU/t = oo, the 
maximum value of the electron concentration is n = 1. 

Now, let us analyze the BCS-BEC crossover in the context of the behaviour 
of the order parameter and the chemical potential. The left panel of Fig. 2 shows 
the dependence of the order parameter at the Fermi level (then, = Aq — {fi/t)Aj) 
on the electron concentration for three values of the intersite interaction. For 
low electron concentrations this is a square root dependence, characteristic of 
the BEC limit. However, for higher electron concentrations one can observe the 
exponential drop of Ak{p,)/t, characteristic of the BCS limit. 

The right panel of Fig. 2 shows the dependence of the modified chemical 
potential on the electron concentration and the intersite interaction (inset). 
As we mentioned before, according to the Leggett criterion [8], the BCS-BEC 
crossover takes place when the modified chemical potential drops below the lower 
band edge (in the 3D case the limiting value is p,/t = —6). The horizontal dashed 
line denotes the band bottom. The values of for n = are exact (despite the 
fact that we have used the Hartree-Fock approximation) and correspond to one 
half of the bound state energy. 

A very important quantity from the point of view of the BCS-BEC crossover 
is the coherence length ^. The suggestions to study this quantity in this con- 
text result from the observations of the properties of exotic superconductors. 
In agreement with the Uemura's plot [10] , the unconventional superconductors 
have a considerably shorter coherence length ^ (~ 20 — 50 A ) than the conven- 
tional superconductors (^ ~ 10^ — 10^ A). The coherence length (pair size) can 
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Figure 3: The dependence of the coherence length in units of the lattice constant 
on the intersite interaction and the electron concentration (inset); the arrows 
show the BCS-BEC crossover points. 



be defined as: 

where <pk = Uki'k = ^k/'^Ek is the pair wave function in the ground state. 

Fig. 3 shows the dependence of the coherence length in units of the lattice 
constant on the intersite interaction and the electron concentration (inset) . Be- 
cause of the fact that we consider the case of hard-core repulsion and intersite 
attraction, ^ exceeds the value of the lattice constant for any interaction values. 
For the values of which slightly exceed the critical value for a given elec- 

tron concentration (5.874 for n = 0), the pair size for small n can be as big as 
the pair size in the BCS regime. When we add electrons, the Pauli Exclusion 
Principle causes a decrease of the pair size, which then becomes just over 1. 
Along with a further increase of the electron concentration (when a crossover 
to the BCS regime for a given interaction strength occurs), the pairs can have 
large sizes {d\ic to the pairing in the momentum space). Fig. 3 shows also 
the dependence of the coherence length on the intersite interaction. For each 
value of n, the asymptote of the function ^ is a line coresponding to the critical 
interaction value for which a transition to the normal state takes place. 

Another quantity which describes the BCS-BEC crossover quantitatively is 
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Figure 4: The dependence of the fraction of condensed particles on the intcrsite 
interaction and electron concentration (inset); the arrows show the BCS-BEC 
crossover points. 



the fraction of condensed particles: 



Fig. 4 shows the dependence of the fraction of condensed particles on the 
intersite interaction and electron concentration (inset). When we increase the 
intersite interaction or decrease the electron concentration, the fraction of con- 
densed particles increases. In the empty lattice limit, no/n takes the largest 
value. The maximum value of the ratio no/n is 0.5, because if all particles in 
the system form pairs and condensate, their concentration is twice lower than 
the initial electron concentration. A very interesting example is the case of 
\W\/t = 5. Because of the fact that this interaction value does not exceed the 
critical interaction value for pair binding, the saturation of the fermion pairs in 
the condensate is very small, as the system is in the BCS regime. 

The case of the fraction of condensed particles exemplifies very well that the 
mean-field approximation fails in the regime of higher electron concentration in 
the hard core repulsion case. Because of the fact that U/t = oo, in the n ^ 1 
limit one should observe a transition to the Mott insulator phase. Indeed, Fig. 
4 shows that the fraction of condensed particles is too high in the n ^ 1 limit. 
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4 Conclusions 



The mechanisms of density-driven and attraction-strength-driven BCS-BEC 
crossover have been considered. The range of occurrence of the local pairs 
phase widens when we decrease the electron conccintration n or we increase the 
attractive interaction \W\. We have considered the s*-wave pairing, for which 
the evolution from the weak coupling limit (BCS) to the limit of the strong 
coupling (EEC) is smooth. For the case of U = 00, both in the 2D and the 3D 
case there exists a critical interaction value for which a bound state forms in 
the empty lattice (when n — > 0) and the critical interaction value for which a 
transition to the normal state takes place. The two critical values coincide in 
the 2D case, but in the 3D case the critical value for the emergence a bound 
state is higher than for the emergence of superconductivity. We have noticed 
that despite the fact that we have used the Hartree-Fock approximation, the 
values of /Z for n = are exact and correspond to one half of the energy of 
the bound state. We have shown that the pair size for small n can be as large 
as the pair size in the BCS regime. Indeed, the Hartree-Fock approximation 
fails for higher electron concentrations in the hard core repulsion case - we have 
considered this issue in the context of the behaviour of the fraction of condensed 
particles. 
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